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ABSTRACT. We characterized the distribution patterns of arboreal lichens in a wet temperate spruce-fir 
(Picea engelmannii/Abies lasiocarpa) rain forest in east-central British Columbia. Using the single rope 
technique for canopy access, we employed the “clump method” to estimate the abundance of arboreal 
lichens in three functional groups: Alectoria, Bryoria, and foliose lichens. Verification was achieved through 
a regression of those estimates with lichen biomass measured on a subset of branches. We used a probability 
proportional to prediction sampling scheme such that the probability that a branch would be selected for 
measurement was proportional to its estimated lichen abundance. Used in this way, the clump method is a 
reliable estimator of lichen abundance that can be used to accurately describe the distribution of lichen 
biomass in tree crowns. 
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RESUMEN. Hemos caracterizado la distribución de muestras de líquenes arbóreos en un abeto-pino (Picea 
engelmannii/Abies lasiocarpa), humedo, templado bosque al este de British Columbia central. Usando el 
método “single rope,” técnica de cuerda simple, aplicamos el método de grupo para estimar la abundancia 
de líquenes abdéreos en tres grupos funcionales; Alectoria, Bryoria y líquenes frondoso. La comprébacion 
fue realizada a través de una regresión de estimaciones con líquenes mediante una prueba aparte (sub- 
sample) obtenido a través de probabilidad proporcional a pruebas pronosticas. La evaluación de la com- 
binacion de esas técnicas muestra que el método grupo (clump method) es un estimador seguro de líquenes 
abundantes que se puede usar para producir un selección datos objectivo que describe con exactitud la 


distribución de epifita líquenes de biomasa. 


INTRODUCTION 


The study of epiphytic lichen biomass in the 
upper canopy of high elevation Engelmann 
spruce-subalpine fir stands has previously been 
limited to studies of felled trees (Edwards et al. 
1960) and studies of lichens within a few meters 
of the forest floor (Stevenson 1979, Van Daele 
and Johnson 1983, Detrick 1985, Rominger et 
al. 1994). This is perhaps due to the difficulties 
involved with accessing the communities in the 
upper canopy coupled with the intimidating task 
of estimating the biomass, which, in some areas, 
exceeds 3000 kg/ha (Edwards et al. 1960). 

Access to the canopy is an integral part of 
studying epiphytic communities. The mass of 
epiphytes on the upper branches is often not ob- 
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vious from the ground (Pike et al. 1975) and an 
inventory of the lower canopy and trunk pro- 
vides an incomplete view of the total ecosystem 
(McCune 1990). Thus, as much effort is often 
dedicated to accessing the canopy as is devoted 
to studying the communities within the canopy 
(Moffett and Lowman 1995). The rewards of 
such effort are found when the canopy organ- 
isms may be studied with as much detail as ter- 
restrial species. This, in turn, increases our con- 
fidence in the reliability of estimation methods 
at each level in the canopy. 

Methods of assessing epiphytic lichen bio- 
mass range from obtaining biomass by tediously 
dissecting and weighing each sample to tradi- 
tional cover estimation techniques using quad- 
rats and transects. Direct biomass sampling is 
too time-consuming to be used in an inventory 
study in which the entire community must be 
quantified. Traditional methods that are gener- 
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ally used in terrestrial plant studies are rapid but 
are not easily applied or verified in a three di- 
mensional environment (Boucher and Stone 
1992, McCune 1990). Quicker, more specialized 
techniques were therefore necessary to describe 
arboreal lichens in this study. We describe the 
use of the “‘clump method,” an estimation meth- 
od that was originally developed to quantify li- 
chen biomass in the lower canopy (Stevenson 
1979), to visually assess lichen abundance 
throughout tree crowns. We also evaluate the ef- 
fectiveness of verifying our estimates with a 
sampling scheme in which the branches with 
high lichen abundance are sampled more inten- 
sively than those with low lichen abundance. 


MATERIALS AND METHODS 
Study Area 


We conducted the research in the Bowron 
River drainage near Pinkerton Mountain in east- 
central British Columbia. The study area 
(53°38'N, 121°25'W) is approximately 90 km 
ESE of Prince George on a gentle south-facing 
slope at 1450—1475 m in the Engelmann Spruce- 
Subalpine Fir Zone (Meidinger and Pojar 1991). 
The climate in the study area is cold and wet. It 
is within the range of woodland caribou (Ran- 
gifer tarandus caribou), a species at risk in east 
central British Columbia that depends on the li- 
chens Bryoria spp. and Alectoria sarmentosa for 
winter forage. 

The coniferous forest is dominated by subal- 
pine fir (Abies lasiocarpa) and Engelmann 
spruce (Picea engelmannii), many of which 
grow in close association with one another, 
forming clumps separated by small natural open- 
ings. Mature trees were 20—30 m tall with char- 
acteristically narrow canopies. Branches were 2— 
4 m long but decrease in length with increasing 
height. They were downwardly sloped and usu- 
ally spaced in clusters of four or more approxi- 
mately every 0.5—1.0 m up the trunk. The den- 
sity of branches increased markedly with in- 
creasing height (Campbell 1998). 


Sample Selection 


The study area is the site of long-term re- 
search on the effects of partial cutting on various 
resource values and ecological processes, in- 
cluding caribou habitat and epiphyte dynamics. 
Study trees were randomly selected from a pool 
of representative trees. To be included in this 
population, trees had to meet the following cri- 
teria: 

a) canopy architecture—trees were located in 
characteristic groups (clumps) of trees, com- 
prising a taller stratum or layer of trees sur- 
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rounded by an unspecified number of shorter 
layers, these clumps physically separated one 
from another; 

b) disturbance vectors—eligible trees were lo- 
cated in undisturbed forest cover, at least two 
tree lengths (ca. 60 m) from man-made open- 
ings (roads, harvest openings); 

c) safety concerns—eligible trees were not in 
immediate proximity to standing dead trees 
(snags) which would pose a hazard to re- 
search personnel and were otherwise deemed 
safe to ascend (no heart rot, major cracks or 
breaks in the main trunk, etc.); 

d) ecological analogues—all trees were located 
within the same subzone of the subalpine 
spruce-fir forest zone, on south-facing slopes 
with less than < 10° inclination, between 
1450 and 1475 m elevation a.s.l, and greater 
than 2 m in height. 


Eight groups of trees were randomly selected 
from within the 4 ha study area that met the 
above noted criteria. A set of 55 trees within 
these clumps made up the sample tree popula- 
tion. Seven additional trees were also included 
in the data set, these representing individual 
trees that had a solitary growth habit (not grow- 
ing in clumps) and were located within 30 m of 
the above selected clumps. All selected trees 
were rigged, climbed, and the amount of lichen 
was estimated. 


Canopy Access 

We reviewed the existing canopy access tech- 
niques (Denison 1973, McCarthy 1988, Perry 
and Williams 1981, Ring and Winchester 1996, 
Tucker and Powell 1991) and judged that most 
were more complicated than was necessary to` 
achieve a safe and reliable working environment 
in our study area. We selected the single rope 
access technique described by Perry (1978) as it 
was the most suitable for the requirements of our 
research in the subalpine canopy where the nar- 
row crowns and dense branches created prob- 
lems in maneuvering around the trunk. 

Using a single rope for ascending and de- 
scending from canopy proved to be quite suc- 
cessful. However, it was not possible to make 
detailed observations in high elevation trees 
from the main climbing line. The branch density 
made it difficult to reach every limb. A second- 
ary system (lanyard) commonly used by tree 
climbers was employed. It allowed us to move 
around to all parts of the tree while remaining 
securely attached to the trunk. 


Epiphyte Biomass Estimation Methods 


Epiphytes were divided into three functional 
groups: Alectoria, Bryoria, and foliose lichens. 
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The Alectoria functional group was composed 
almost entirely of Alectoria sarmentosa, al- 
though Usnea spp., which are extremely rare in 
the study area, were included if they were pre- 
sent. The Bryoria functional group comprised 
Nodobryoria oregana and a number of Bryoria 
species, primarily Bryoria capillaris, B. fremon- 
tii, B. fuscescens, B. glabra, and B. pseudofus- 
cescens. The foliose lichen functional group was 
dominated by Hypogymnia imshaugii, H. me- 
taphysodes, H. occidentalis, H. physodes, H. tu- 
bulosa, Parmelia sulcata, Platismatia glauca, 
and Tuckermannopsis platyphylla. 

Those groups were selected because of their 
functional and ecological importance in the 
study area. Caribou use both Alectoria and Bryo- 
ria, but prefer the latter (Rominger et al. 1996). 
The two fruticose groups differ in their distri- 
bution in the canopy, and are expected to re- 
spond differently to forestry activities that alter 
the structure of the canopy. Foliose lichens con- 
taining green algal symbionts are abundant in 
the study area, but are not known as important 
forage items for wildlife. Epiphytic bryophytes 
and cyanolichens were not included because 
they are rare or absent in our study area. Crus- 
tose lichens were excluded because of time con- 
straints. | 

The abundance of each functional group on 
each branch was estimated using the “‘clump 
method,” a technique that facilitates estimation 
of biomass on a three dimensional substrate 
(Stevenson 1979). The lichen on a given branch 
was compared to a standard unit of lichen of 
known size and weight; biomass was estimated 
according to how many multiples of the standard 
unit that it represented (FIGURE 1). For both fru- 
ticose lichen functional groups (Alectoria and 
Bryoria) a 2.5 g clump of lichen was used as a 
standard and the total biomass of fruticose li- 
chen was estimated as multiples of that weight. 
Alectoria and Bryoria were subsequently parti- 
tioned into separate functional groups by giving 
an estimate of the percentage of the total fruti- 
cose lichen that was made up of Bryoria. A sim- 
ilar method was used for the foliose lichen func- 
tional group. The three dimensional clump of 
fruticose lichen was modified to a “‘card”’ of fo- 
liose lichen of known dimensions. The standard 
was an 8.5 X 8 cm surface which was the equiv- 
alent of approximately 1.5 g of foliose lichen. 
The foliose lichen on each branch was then as- 
sessed according to how many multiples of that 
surface area (card) it represented. In many cases 
foliose lichen grew in horizontal layers, forming 
mats. In such cases, the number of cards re- 
quired to make the mat was estimated by the 
observer. 

All biomass estimations in one tree were 








The standard clump and card used for 
lichen biomass estimation at Pinkerton Mountain. The 
clump of fruticose lichen (Alectoria and Bryoria) is 
2.5 g, of which 30% is Bryoria and the card of foliose 
lichen weighs 1.5 g. 


FIGURE 1. 


made by the same researcher to minimize inter- 
observer error and to maintain the high level of 
consistency required with this technique. As 
well, observers calibrated their estimates 
throughout the field season by practicing togeth- 
er. Accuracy and precision were tested by com- 
paring the estimates to dissected samples ob- 
tained by sampling with a probability of selec- 
tion that was proportional to estimated biomass. 


Verification of Estimates 


Because biomass sampling is time-consum- 
ing, we aimed at a sample size of 5% of the 
branches on each sampled tree. Random sam- 
pling with a relatively small sample size is un- 
likely to adequately represent the branches with 
high lichen biomass, because the frequency dis- 
tribution of lichen biomass on branches is 
strongly skewed right. We used Probability Pro- 
portional to Prediction (3P) sampling (Cochrane 
1977), a method similar to Probability Propor- 
tional to Size (PPS) sampling as described by 
Pike et al. This provided better representation of 
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branches with high biomass loading, as com- 
pared to a strictly random sampling of canopy 
branches. Branches were chosen for inclusion in 
the sample in direct proportion to the estimated 
biomass of lichen found on them. The resulting 
sample thus has a greater proportion of branches 
with a large amount of lichen than would occur 
through random sampling. 

Branch selection was achieved in three steps. 
First, the number of clumps of fruticose lichen 
on each branch of a tree was estimated, as de- 
scribed above. Second, the sample branches for 
that tree were selected, using a method borrowed 
from timber cruising (Dilworth and Bell 1971). 
For each branch, a random number was drawn 
between 1 and an upper limit. The formula for 
choosing that top random number, T, was: 


T = Sx/s 


where T was the top random number, Sx was the 
sum of the estimates of total fruticose lichen for 
the entire tree, and s was the expected sample 
size, which was 5% of the branches on the tree. 
The estimate of fruticose lichen on each branch 
was then compared to a random number less 
than or equal to the top random number. The 
branch was selected for sampling if the estimat- 
ed biomass was greater than or equal to the ran- 
dom number. 

Sometimes, because of the skewed frequency 
distribution of the lichen estimates, T was small- 
er than one or more of the largest single esti- 
mates on a tree. In that case, the formula was 
modified so that T was doubled and each branch 
on that tree was given two chances to be sam- 
pled. If a branch was selected by both random 
numbers then it was counted twice in the cal- 
culations (Stevenson 1979). The third step was 
to return to the tree and remove all the lichen 
from the sample branches. 

Natural log (Bryoria +1) and log (Alectoria 
+1 and foliose +1) transformations were ap- 
plied to the measured and estimated lichen abun- 
dance values, to which a constant of 1 had been 
added to avoid taking the logarithm of zero. 
Simple regression models were calculated for 
each functional group to describe the relation- 
ship between the measured and estimated lichen 
on the 5% sample. The resulting (In-In or log- 
log) regression equations were used to calculate 
a corrected lichen biomass estimate for each 
functional group on each branch on each study 
tree, including the trees that had not been sub- 
jected to destructive sampling. Estimates were 
converted back to g/branch by taking the anti 
log/In of the output of the regression equation. 
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RESULTS 


Our measurements of lichen biomass were 
based on a combination of the clump method 
and 3P sampling. The standard errors of estimate 
(SEM) indicate that there was some variability 
in the difference between estimated and mea- 
sured biomass. However, the strong coefficients 
of determination for the regressions between es- 
timated and measured biomass (FIGURE 2) and 
large effect sizes for each functional group 
(Alectoria: R? = 0.68, f? = 2.12, SEM = 6.41; 
Bryoria: R? = 0.76, f? = 3.14, SEM = 19.63; 
foliose lichen: R = 0.61, f? = 1.58, SEM = 
7.53) showed that the clump method coupled 
with 3P sampling for verification was an effec- 
tive tool for assessing lichen biomass at the tree 
level. 

There were no significant differences in the 
quality of the estimates on branches measured 
early versus late in the season or on branches 
with low versus high abundance (R? difference 
test, P < 0.05; Fox 1997). That indicates that 
the clump method can be effectively applied 
with little practice time throughout the canopy. 
It was, however, necessary to “‘calibrate”’ the es- 
timates of the observers at the beginning of the 
field season, based on side-by-side comparisons 
using standardized lichen clumps of known 
weight, to ensure consistency within the re- 
search team. 

The application of the clump estimation tech- 
niques in sub-alpine spruce-fir forest at Pinker- 
ton Mountain demonstrates a strong height re- 
lated vertical zonation for each of the three func- 
tion groups (FIGURE 3). Members of the Bryoria 
functional group reach greatest abundance in the 
upper canopy, up to 30 g/branch, at heights 
above 15 m in the canopy. Alectoria thalli, in 
contrast show greatest abundance in the lower 
canopy, ca. 5 m above the forest floor. Foliose 
lichens are represented are represented at all 
height stratums, although reaching greatest 
abundance in the lower canopy. 


DISCUSSION 


The single rope technique was successfully 
applied in the high elevation forests of east cen- 
tral British Columbia. The equipment was rela- 
tively inexpensive and readily available from 
many outdoor equipment stores, could be carried 
and easily assembled by one person, and afford- 
ed rapid access. 

The success of the access techniques allowed 
us to apply the clump method confidently to 
each branch individually. Stevenson (1979) ap- 
plied the clump method to both highly visible 
and poorly visible branches and found that co- 
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FIGURE 2. Regression lines for the clump method 
as an estimator of lichen abundance. Measured lichen 
(Y) is the biomass of the dissected and weighed sam- 
ples obtained through 3P sampling, and estimated li- 
chen (X) is the abundance of lichen estimated using 
the clump method on the same sample as above. Re- 
gression statistics were performed on log and In trans- 
formed data. A. Alectoria: Log Y = 0.0755 + 0.992 
log X (R? (N = 94) = 67.9%). B. Bryoria: Ln Y = 
0.790 + 1.102 In X (R (N = 94) = 76.1%). (C) Fo- 
liose: Log Y = 0.136 + 0.969 log X (R? (N = 94) = 
61.3%). 


efficients of variation ranged from 15% to 85%. 
Coefficients were lowest for branches that were 
most easily observed and not concealed by oth- 
ers. Combining the clump method with tree 
climbing techniques ensured that every branch 
in this study was highly visible. 

The single rope technique provided near total 
access to the tree crown. As with other canopy 
studies (Pike et al. 1977), safety concerns pre- 
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FIGURE 3. Height related vertical zonation of li- 
chen biomass at 1 m class intervals (Mean and Stan- 
dard Error) for Alectoria, Bryoria, and Foliose func- 
tional groups in subalpine spruce-fir forest at Pinkerton 


Mountain, B.C. 


vented the evaluation of the very top of each 
tree. However, since the branches in this region 
were primarily young, most had a diameter be- 
low our study cut off (2 cm) and would not have 
been included in the sample. Thus, this omission 
only slightly affected our overall lichen biomass 
estimates. In addition, since the vast majority of 
lichen in that portion of the tree was Bryoria, 
the impact of excluding these branches on Alec- 
toria and foliose lichen was negligible. 

We found that the clump method was well 
suited for the lichen community structure of our 
study site in two ways. It allowed rapid esti- 
mations of all three lichen functional groups on 
each branch. This was necessary both because 
this was an inventory study and because our 
scheme for verifying and correcting our esti- 
mates required that the entire population be as- 
sessed. Also, we found that even though esti- 
mates could vary by as much as two orders of 
magnitude, we were able to use a single refer- 
ence clump for fruticose lichens. We were also 
able to use a single card as a standard for foliose 
lichens in this ecosystem in which the branches 
were small and easily seen. Although our eco- 
system required only one size standard, a refer- 
ence clump can be tailored in size to suit the 
community, and several reference clumps of dif- 
ferent sizes can be used if necessary. This flex- 
ibility makes the method applicable to a wide 
range of ecosystems. 

The efficacy of the clump method as an esti- 
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mator varied among our lichen functional 
groups. There are two reasons why this variation 
might exist. First, the standard clump was large 
and therefore more easily compared to lichen in 
large amounts. Bryoria makes up the bulk of 
epiphyte abundance and is therefore more ac- 
curately estimated using a large standard than 
either of the less abundant functional groups. 
Second, the 3P sampling was calculated based 
on the estimates of total fruticose lichen (Alec- 
toria and Bryoria combined), the majority of 
which is Bryoria. In essence, then, the selection 
of samples was based on the sum of estimates 
of Bryoria and therefore has a larger proportion 
of branches heavily laden with Bryoria than 
those with large loads of either Alectoria or fo- 
liose lichen. The larger coefficient of determi- 
nation for Bryoria reflects this bias. Sampling 
according to Bryoria abundance was, however, 
still a reasonable estimator of both Alectoria (R? 
= 67.9%) and foliose (R? = 61.3%). It appears 
advantageous to select branches for sampling ac- 
cording to the most important or most abundant 
attribute in the context of the study (i.e., caribou 
forage lichen). 

By estimating lichen abundance on all branch- 
es, measuring lichen abundance on some 
branches, and using the relationship between es- 
timated and measured values to correct the es- 
timates, we arrived at a reasonably accurate and 
objectively collected data set. The variation in 
the relationship between estimated and mea- 
sured biomass, which was reflected in the SEM, 
is a weakness in any epiphytic lichen estimation 
technique that occurs because it is not logisti- 
cally possible to directly measure every branch 
in the population. Our technique resulted in a 
relationship between measured and estimated li- 
chen biomass that was comparable to those re- 
sulting from many of the existing methods. The 
coefficients of determination from the regression 
equations calculated for the three functional 
groups varied from 0.76 for Bryoria to 0.61 for 
foliose lichens. McCune (1990) estimated lichen 
biomass using a variety of estimators. Resulting 
coefficients of determination were 0.58 < R < 
0.91, 0.30 < R? < 0.96 and 0.33 < R? < 0.91 
using cover classes, maximum thallus length, 
and maximum thallus width, respectively. Al- 
though our estimator was not as closely related 
to the actual lichen biomass as in the best cases 
with previously used estimators, it is uncertain 
whether cover classes and maximum measures 
for each functional group would be easily ap- 
plied in ecosystems where Alectoria and Bryoria 
thalli are abundant and intertwined (Stevenson 
and Enns 1993). 

Regressions, rather than ratios (Pike 1981), 
were used in this study to compare the estimates 
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with the sub-sampled biomass. The ratio method 
is more suitable for intensive studies of one or 
a few trees, because correction ratios are calcu- 
lated for the whole population that is sampled 
with probability proportional to prediction. In 
studies that involve climbing, that population 
would normally be a single tree. It would be 
impractical to estimate abundance on the entire 
population of trees, select the sample branches, 
and return to the trees to collect the samples. 
Thus, the number of branches sampled on each 
individual tree would have to be sufficient to 
achieve a reasonable confidence interval when 
using ratios. By using regression statistics, we 
were able to pool the data from all our study 
trees, and apply a single regression for each 
functional group to all branches. 

The efficacy of these clump estimation tech- 
niques is demonstrated by the fine resolution 
data obtained for height related vertical zonation 
of lichen functional groups in subalpine forest 
at Pinkerton Mountain. Lichens in the Bryoria 
functional group reach their greatest abundance 
in the upper canopy, in marked contrast to that 
of the Alectoria functional group lichens. Al- 
though Alectoria and Bryoria spp. lichens have 
previously been grouped together within an 
Alectorioid functional group (e.g., McCune 
1990) our present findings demonstrate quite dif- 
ferent ecological affinities for these two genera, 
at least within the context of these moist-cool 
subalpine forests. 

This differential response of these two lichen 
groups may reflect changes in both substrate 
characteristics and operating environments for 
lichen growth, especially parameters of canopy 
microclimate. Goward (1998) hypothesizes that 
Bryoria lichen communities cannot tolerate ex- 
tended wetting episodes, and thus are favored in 
upper canopy aspects. We are presently profiling 
canopy microclimate in context of measures of 
lichen physiology to verify the applicability of 
this hypothesis. A more detailed examination of 
canopy structural variables for this site is pro- 
vided in Campbell (1998). 

In summary, the combination of the clump 
method, 3P sampling, and regression statistics 
allowed us to characterize lichen abundance and 
distribution on a large number of trees. These 
techniques may be effectively applied where de- 
scriptions of lichen biomass at the stand level 
are required. In the subalpine spruce-fir forests 
of the present study area these measures show 
that Alectoria and Bryoria lichens occupy dif- 
ferent niches within the canopy, and suggest that 
traditional groupings of these genera into a com- 
mon Alectoriod functional group may not be ap- 
propriate in these cool wet ecosystems. 
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